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The intermolecular force constants for a 6-12 Lennard-Jones potential model have been 
computed for air from second virial coefficients derived from experimental pressure- volume- 
temperature measurements. Nitrogen: oxygen interaction second virial coefficients and 
interaction force constants have also been derived. 



The following symbols are used throughout the 
paper: 

Z = Compressibility factor, Pv/RT, dimensionless. 
P= Absolute pressure. 
v= Molar volume, cm 3 /mole. 
R = Universal gas constant, in units of Pv/T. 
T= Absolute temperature. The triple point of water is 

defined as T=273.16°K. 
B = Second virial coefficient, cmVmole. 
C= Third virial coefficient, cm 6 /mole 2 . 
D = Fourth virial coefficient, cm 9 /niolc 3^ 
E= Fifth virial coefficient, cm 12 /mole 4 . 
€ = Maximum energy of binding between a pair of molecules 

with a Lennard-Jones potential. 
k= Boltzmann's constant. 
AT t =Mole fraction of the ith component. 
b= Lennard-Jones interaction parameter, § Nird 3 , in units of 

molar volume. 
N = Avagadro's number. 
d= Collision diameter of a pair of molecules at zero energy 

of Lennard-Jones interaction. 
A= Angstrom units (=cm~ 8 ). 
e= Potential energy of two interacting molecules. 
r — Separation between centers of two interacting molecules. 
r = Collision diameter of a pair of molecules at maximum 

energy, e , of Lennard-Jones interaction. 

Compressibility factors for dry air have been 
determined experimentally and are reported else- 
where [l]. 2 By fitting the P-v-T data to a Kamer- 
lingh-Onnes virial equation of state 



~'RT 



= 1 + B/v+C/v 2 +D/v*+E/v* 



(1) 



the virial coefficients B, C, D, and E were determined 
at various temperatures in the region between the 
critical and 0° C. These are tabulated in table 1. 

The second virial coefficients of nonpolar mixtures 
can be derived, to a good approximation, from the 
force constants of the pure components [2, 3]. If we 
assume that the intermolecular forces of the pure 
constituents of dry air (N" 2 , 2 , and A) can be de- 
scribed by a Lennard-Jones model, then we may 
define pseudoforce constants for the air mixture: 



(€o/k) m i x =^2N t (€olk) t 



(2) 
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and 



W=§ *#(<*)£*=! tN^NM 3 . (3) 



If we assume that the air mixture has the com- 
position 78.2 percent of N 2 , 20.8 percent of 2 , and 
1.0 percent of A (cf. [1]), and the N" 2 , 2 , and A can 
be represented by a 6-12 Lennard-Jones intermolecu- 
lar potential function with force constants: 



and 



(e /k) N =95M o , 
(€o/&)o 2 =118.0°, 

(€ /&)a=119.5°, 



(rf) N =3.69 A [3, 4] 



(d)o a =3.46A [5] 



(d) A =3.42 A, [6] 



then, from eq (2) and (3), (e /k) &lr = 100.8° and 6 alr = 
G0.82 cm 3 /mole. 

These 6-12 Lennard-Jones pseudoforce constants 
can then be used to compute 4 the values of the second 
virial coefficients for air at various temperatures. 
The values thus computed are in good agreement 
with the experimentally determined second virial 
coefficients [1]. 

Table 2 compares the calculated and experimental 
values of B. 

The second virial coefficient of a binary gas mix- 
ture can be expressed in terms of the second virial 
coefficients of the constituent gases and an inter- 
action coefficient related to collisions of the unlike 
molecules. 

BmIx=SSAWB„. (4) 

i j 

We may then write for air, noting that B n =B 2 \ ) 

B &lT = (0.79) 2 7? N2 + (0.21) 2 £o 2 +0.33187A 2 . (5) 

Rearranging: 

B 12 =B N2 : = (7? air -0.624l7^ 2 -0.044l5o 2 )3.014 (6) 

Using the P-v-T data of Friedman et aL, for B alT 
[1], of Friedman for S Na [3], and the tabulation of 
Hilsenrath et al., for Bq 2 [6], the N 2 :0 2 interaction 
second virial coefficients as a function of tempera- 
ture may be calculated from eq (6). These arc 
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Table 1. Virial coefficients of air 



T 


£* 


C 


D 


E 


° K. 


cm 3 lmole 


cm 6 /mcle 2 


cm^/mole^ 


cm^/moleA 


273. 15 


-12.57 


0.7343X103 


0. 1468X10« 


-0.0688X10 8 


250 


-19.30 


.7646 


.1836 


-.1011 


200 


-38. 20 


1.6303 


-.0116 


.0345 


175 


-55.54 


2. 5492 


-.0442 


.01986 


150 


-76.36 


3. 3923 


-.1101 


. 04184 



a The estimated standard deviation of the second virial coefficients, B, vary 
from 4 to 15 percent of the value of the coefficients. The standard deviations of the 
other coefficients are a large fraction of the value of the coefficients. 

Table 2. Second virial coefficient of air 



T 


B 








Calculated 


Experi- 
mental 


K 


cm 3 lmole 


cmz/mole 


273. 15 


-13.4 


-12.57 


250 


-19.6 


-19.30 


200 


-38.9 


-38. 20 


175 


-53.3 


-55. 54 


150 


-74.1 


-76.36 



Table 3. Nitrogen: oxygen 
interaction second virial 
coefficients 



T 


-Bn 2 :o 2 


K 


cm^/mole 


273. 15 


-17.5 


250 


-24.0 


200 


-44.0 


175 


-61.3 


150 


-84.7 



tabulated in table 3. These interaction virial 
coefficients can be fitted to a Lennard- Jones 6:12 
potential model. The intermolecular force constants 
for the interaction between a nitrogen molecule and 
an oxygen molecule thus computed are: 

eo/k= 106.4° 

b= 61.1 cm 3 /mole 

d= 3.64 A. 

The intermolecular potential energy is plotted as a 
function of the separation of the centers of the nitro- 




Inter molecular potential between a nitrogen mole- 
cule and an oxygen molecule (solid line). 

The potential energy curves for oxygen-oxygen interactions and for nitrogen- 
nitrogen interactions are indicated by the dashed lines. 

gen from the ox} T gen molecules in figure 1. The 
equation for this curve is 

H(f)[(F-(f] 



or 



^=425.6[(3.64/r) 12 -(3.64/r) 6 ]. 



(8) 



The collision diameter at the maximum energy of 
binding between the pair of interacting molecules 

is r o = d^/2 . For nitrogen : oxygen collisions, r = 4.08 A. 
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